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Abstract. The Medusa algorithm takes as input two postcritically fi- 
nite quadratic polynomials and outputs the quadratic rational map 
^ which is the mating of the two polynomials (if it exists). Specifically, the 

output is a sequence of approximations for the parameters of the rational 
map, as well as an image of its Julia set. Whether these approximations 
converge is answered using Thurston's topological characterization of 
rational maps. 

This algorithm was designed by John Hamal Hubbard, and imple- 
mented in 1998 by Christian Henriksen and REU students David Farris, 
j^/y and Kuon Ju Liu. In this paper we describe the algorithm and its im- 

plementation, discuss some output from the program (including many 
• pictures) and related questions. Specifically, we include images and a 

f^h discussion for some shared matings, Lattes examples, and tuning se- 

^ quences of matings. 

B 

1. Introduction 

^ The study of the dynamics of rational maps of the Rieman sphere is 

^ greatly facilitated by the fact that a wide variety of dynamical phenomena 

can be illustrated using only the quadratic family P c (z) = z 2 + c. Of course 
(N| most general theorems about rational maps have examples in the quadratic 

family, but further, in some cases the dynamics of a quadratic polynomial 
appear within a rational map. The most basic example of this phenomena 
• • is through polynomial-like behavior. In addition, there are several ways 

to combine two (or more) quadratic polynomials to produce rational maps 
^ whose dynamics can be described via a combination of the quadratic poly- 

nomial dynamics. Probably the first such example was a polynomial mating 
discovered by Adrien Douady |Dou83| . 

In order to define matings, first we must step back to quadratic poly- 
nomials. It is simple to write a computer program which, given a c, will 
compute (approximately) the orbit of any given point under the quadratic 
polynomial P c . To illustrate the overall behavior one draws the filled Julia 
set, K C1 the set of points whose orbit under P c does not tend to oc. This 
also illustrates the Julia set, J c , the topological boundary of K. (See 
Figure [I] for a sample J c .) We may examine experimentally the dynamics of 
one map at a time with such a program. 
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The next natural step is to understand how the dynamics changes with 
a change in the parameter, c. We organize the parameter space by defining 
M, the Mandelbrot set, as the set of all c in C for which the Julia set J c is 
connected (see Figure [3]). By Fatou's fundamental dichotomy theorem, 
this is equivalent to the set of all c such that the orbit of the critical point 
under P c lies in K c . Thus it is also a simple matter to generate a picture 
of M, and a program which will draw the Julia set J c when a parameter c 
in M is selected. After a brief investigation with such a program, one sees 
intriguing patterns, and a relationship between M and the Julia sets of its 
children, the quadratic polynomials. 

In addition to the definition of M, many basic results in the theory of the 
iteration of rational functions support the premise that the behavior of the 
critical orbit is crucial for describing the dynamics. The dynamics are most 
amenable to analysis when the polynomial P c is post critically finite (PCF), 
i.e., the orbit of the critical point is finite. A key technique in giving 
a mathematical description of the patterns of quadratic polynomials turns 
out to be combinatorics. For a postcritically finite quadratic polynomial, 
we can build a labelled graph, called a spider, which gives a combinatorial 



description of the dynamics of the polynomial. This is described in {2.2 



The reverse problem, of starting with a combinatorial spider and pro- 
ducing a quadratic polynomial P c (i.e., producing a parameter c) whose 
dynamics are given by that model, is solved by the spider algorithm. The 
spider algorithm is an iterative procedure, based on Thurston's topological 
characterization of rational maps |DH93j . and is described fully in [HS94j. 

The main subject of this paper is the Medusa algorithm, which takes 
two combinatorial spiders, glues them together in a certain manner (hence 
the name Medusa), then runs a sort of double spider algorithm which, if it 
converges, produces a rational map which is the mating of the two quadratic 
polynomials associated with the originally inputted spiders, see Theorem |3.9| 

John Hamal Hubbard designed the Medusa algorithm, based on Thurston's 
theory QDH93J) and the foundational theory of polynomial matings devel- 
oped by Douady, Hubbard, Shishikura, Rees, Tan Lei and others ( |Dou83t 
Ree92l ILei921 IShiOOj . see SlOb. The computer program implementing the 
algorithm was written under Hubbard's direction by David Farris, Chris- 
tian Henriksen and Kuon Ju Liu, in a 1998 summer research experience for 
undergraduates program. The full source code for Medusa is available for 
download at |Dyn |. 

Some progress has been made in the study of polynomial matings since 
1998, however there are still many intriguing questions. The goals of experi- 
mental software like Medusa are to help form conjectural answers to existing 
questions, as well as inspire new questions. After explaining the algorithm 
and implementation, in the final section of this paper we provide several 
examples of images we created using Medusa, which serve to illustrate and 
examine several of the phenomena of matings. Specifically, we include im- 
ages and a discussion for some Lattes examples, shared matings, and tuning 
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sequences of matings. We hope this paper will energize future researchers to 
study polynomial matings, and we expect Medusa is of service in advancing 
the field. 

Organization of sections. In ^2] we provide needed prerequisite material 
on the dynamics of quadratic polynomials and polynomial matings. In 
we describe the Medusa algorithm and its implementation, proving Theo- 
rem 3.9 The final section, £|4j contains examples of output from the program 
related to a few areas of interest in the study of matings. 

Acknowledgements. The authors thank Dierk Schleicher, Adam Epstien 
and Tan Lei for inspiring discussions and advice on how to write this paper. 
All images of Julia sets of quadratic polynomials were generated with the 
Otis fractal program |Kawj . 

2. Background 

2.1. Notation. We write C = C U {oc} for the Riemann sphere, i.e., the 
one point compactification of the complex plane, endowed with the complex 
structure with respect to which the identity restricted to C is a chart, and 
z i y \ I z <\> conformal isomorphism. We write S 2 for C viewed as a topological 
manifold, i.e., not equipped with a canonical complex structure. 

2.2. Quadratic polynomials and combinatorics. If K c is connected, 
then there is a unique conformal isomorphism 

tp c : C-D^C-K C , 

such that ^c(oo) = 1. This map conjugates w \-> w d to P c . The curve 
TZt(c) — TZt — {^(re 27Tlt ) : r > 1} is the external ray of angle t. For a 
postcritically finite polynomial the filled in Julia set K c is locally connected, 
and then ^ c extends continuously to the boundary. If we parameterize the 
circle by R/Z, then the map tp c on the boundary becomes 

7c :R/Z^ J c , 

and 7 C is a semiconjugacy of multiplication by two to /, i.e., 7 c (2t) = 
Pcilit)). Then 7 c (t) is called the landing point of TZ t {c). Call 7 C the 
Caratheodory map of P c . See Figure [I] for a picture of a Julia set and 
some external rays. 

Given a postcritically finite quadratic polynomial, P c , choose 9 C E R/Z so 
that 1Zq c is the external ray associated with the critical value, c. That is, 1Zq c 
lands at c, if c E J c - Otherwise the critical point is periodic. If the critical 
point is fixed, take 6 C = 0. If the critical point is periodic of period n > 1, 
the critical value is contained in the immediate basin U of a superattracting 
cycle and there exists a pair of rays landing at the root of U whose closure 
seperates the critical value from the other points in the critical orbit. Take 
6 C to be one of the two angles corresponding to this pair of rays. 

We can use 7 C and 9 C to create a simple combinatorial model of the critical 
orbit. 
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Figure 1. Left: the Julia set of z H> z 2 + z, and critical 
orbit rays 1/6, 1/3,2/3. Right: the Julia set of z H> z 2 — 1, 
and critical orbit rays 1/3,2/3, plus 1/6 for comparison. 



Given a rational number 9 E R/Z, following Hubbard and Schleicher 
QHS94J) we define the standard 9 -spider, S^cC by: 

§, = {re 2 - 23 " l0 :r>l,j = l,2,...}U{oo}. 

See the image on the left in Figure [2] for an example, it shows the spider for 
one of the Julia sets of Figure [Tj One may view this as a spider, with legs 
the rays emanating from the unit circle which are in the orbit of 9 under 
angle doubling, and body the point at infinity. 

Since 7 C semi- conjugates P c to angle doubling, 7 C maps S>q c to the union 
of 1Zq c and its images under P c , plus the point at infinity. Note if 9 is 
rational, then it has finite orbit under angle doubling, so the spider has a 
finite number of legs. Similarly, if P c is postcritically finite, then 9 C will be 
rational. We denote the endpoints on the unit circle of the spider legs by 
Z j = e 2Mj ~ le . 

The spider illustrates the critical orbit. Using this diagram we can also 
create a sequence called the kneading sequence of 9 which records informa- 
tion about the order of the critical orbit in this diagram. Take the plane 
containing the spider and cut along the line composed by the rays of 
angle 9/2 and (9 + l)/2. Label by A the open half of the plane containing 
9, label the other open half B. See the right hand image of Figure |2| Label 
the ray of angle 9 by * a , and the ray of angle (9 + l)/2 by For any angle 
t, its ^-itinerary is the infinite sequence of labels from (A, B, * a , *&) corre- 
sponding to the position in the labelled plane of the points in the forward 
orbit of t under angle-doubling. The kneading sequence of 9, denoted k(9), 
is the ^-itinerary of the angle 9. Note a symbol * n appears in this sequence 
if and only if 9 is periodic under angle doubling. 
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Figure 2. Left: the spider for 9 = 1/6. The critical orbit 
is (1/6 i — ^ 1/3h>2/3h-> 1/3). Right: the kneading sequence 
for this spider is if (1/6) = A AB. This spider models f(z) = 
z 2 + z, whose Julia set is shown in Figure [TJ 

In this paper, we are interested in combining and comparing quadratic 
polynomials. In order to keep track of the dynamics of various maps we 
are studying, we use the discovery of Douady and Hubbard (see [DH82]) on 
how 9 C relates to the position of c in the Mandelbrot set. They show the 
Mandelbrot set, M, is connected, with simply connected complement in C, 
hence there is a unique conformal isomorphism ^ m • C — M — >► C — D which 
fixes oc and such that ^^(oc) = 1. Then ^ m defines external rays outside 
of M, by images of straight rays outside of the disk. It happens that for 
any rational angle 9 — p/q, the map $m extends radially to the boundary, 
to define a landing point c{9) for the ray of angle 9. Given a postcritically 
finite polynomial P c to which we associate the angle 9 C , then the parameter 
ray of angle 9 C will either land at c (in the preperiodic case) or at the root of 
the hyperbolic component of M that has c as a center (in the periodic case). 
For example, for the basilica, f(z) = z 2 — 1, the external rays associated 
with the critical value —1 is of angle 1/3 and 2/3. The parameter rays of 
angle 1/3 and 2/3 lands on the Mandelbrot set at the root point of the bulb 
containing the basillica (the real bulb). Figure [3] shows the Mandelbrot set 
and some external rays. 

2.3. Mating quadratic polynomials. Let f n (z) = z 2 + c n , n — 1, 2 be 

two quadratic polynomials, with Julia sets J n . Assume each J n is locally 
connected, and 7 n is the Caratheodory map of f n . Define /C = K\ U K2/ ^ 
to be the quotient space of the disjoint union of K\ and K2 in which for 
each t E R/Z, we identify 71 (t) with 72 (— t). In other words, we obtain a 
topological space K by gluing K\ and K2 together along their boundaries 
via 71 (t) ^ 72 ( — t). Consider this definition while viewing Figure [!} In 
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Figure 3. The Mandelbrot set, i.e., the set of all 
c in C for which the Julia set J c is connected, 
shown above in black, together with the external rays: 
0, 1/511, 1/7, 10/63, 1/6, 3/14, 1/5, 1/4, 169/511, 1/3, 
255/511,1/2,2/3,5/6. 

general one might imagine K as some bizarre balloon animal (possibly with 
infinitely many body segments), but we will see below that in many cases, 
/C is simply a sphere. 

On the space /C, define the map f\ JL / 2 by f n on K n , n = 1, 2. Since 7 n 
semiconjugates / to multiplication by two on J n , this map is well-defined 
and continuous (no matter how bizarre the space K may be). 

If there is a quadratic rational map F which is topologically conjugate on 
C to fi JL / 2 on /C, then F is called a mating of /i and / 2 . We denote this 
relationship by F = fi JL / 2 , and in this case say the mating of /i and / 2 
exists. The conjugacy /i: /C — > C is required to be an orientation preserving 
homeomorphism which is holomorphic on the interiors of each K n . It is 
believed that if F exists, it is unique up to Mobius conjugation. 

Note that a mating of any quadratic polynomial f\ with f2(z) — z 2 yields 
F = fv 

Results of Rees, Shishikura, and Tan Lei ([Ree92 ; ILei92[ IShiOO] ) show 
that whether the mating of two PCF quadratic polynomials fi and / 2 exists 
can be answered in terms of the location of c\ and c 2 in parameter space. 
The fundamental existence theorem is: 

Theorem 2.1. are PCF quadratic polynomials, TFAE: 

• K is homeomorphic to the sphere S 2 ; 
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• there exists a quadratic rational map F which is the mating of f\ 
and f 2 ; 

• c\ and C2 do not belong to complex conjugate limbs of the Mandelbrot 
set, M. 

We refer the reader to Milnor's book [Mil99j for detailed background 
on the dynamics of polynomial maps of C, and his article [Mil04] for a 
more complete discussion of the definition of mating and its subtleties, a 
discussion of many foundational results on matings, and a detailed analysis 
of an interesting example of mating. 

3. From Thurston's Algorithm to the Medusa Algorithm 

Thurston's algorithm is a proof that given a branched covering g of the 
sphere there exists a rational map F that is Thurston equivalent to g unless 
there exists a Thurston obstruction. The proof can be made into an iterative 
procedure computing a sequence of complex structures and rational maps F n 
which, when properly normalized, converges to F. In this section we see that 
we can take g to be a model of the mating of two quadratic rational maps, 
and extract finite dimensional but crucial information about the complex 
structures produced by Thurston's Algorithm so that the sequence F n can 
be recovered. This is the heart of the Medusa Algorithm. Because of the 
finite dimensional information needed to run the algorithm, it lends itself to 
actual computation. 

3.1. The Theory. 

Normalizing matings. Assume /i,/2 are postcritically finite quadratic 
polynomials and F = f\ _u_ f 2 . Each f n has one critical point 0, which lies 
in K n . Thus F has two distinct critical points. By conjugating F with a 
mobius transformation we can arrange that the critical point coming from f\ 
is at the origin, the other critical point at infinity and the two glued-together 
beta fixed points are at 1. Therefore we know that any such mating belongs 
to the following family of maps. 

Notation 3.1. We normalize the rational maps which are matings by: 
(1) T — {F rational of degree 2 | 0, oc are critical points and F(l) = 1}. 

Note that every rational map of degree two is conjugate to (at least one) 
member of T . 

The following innocent lemma, which is trivial to prove, is of fundamental 
importance to why there is such a thing as the Medusa Algorithm. 

Lemma 3.2. Given two distinct points u, v G C \ {1} there exists a unique 
F E T so that F(0) = u and F{o6) = v. 

The lemma shows that there is some magic to quadratic rational maps. 
Normalized in the way described, we just need the position of the two critical 
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values (and which correspond to which critical point) to uniquely determine 
the map. We don't need any extra combinatorial information. 

Proof. We prove the lemma in the case where u, v are different from infinity. 
The case where either u or v equals infinity is just as easy and left to the 
reader. First notice that F : z ^ ^ u ~ 1 \ vz 2 ~j( v ~^ e J 7 has the desired 

(u— l)z z — [v— 1) ' 

properties, so we need to show that this is the only such map in T . Since 
the origin and infinity are critical points, we can write 



F(z) 



az 2 + c 
bz 2 + d' 



That 1 is fixed, F(oo) — v and F(0) — u implies that a — vb = 0, c — ud — 
and a — b + c — d = 0. When either u or v is different from 1 the matrix 

1 -v 
1 -u 
1-11-1 

has rank 3. It follows that every solution to the three equations can be 
written (a, 6, c, d) — \{{u — l)v, u — 1, — u(v — 1), —(v — 1)) for some A E C, 
and therefore F is uniquely determined. □ 

In the following we will write F u v for the map given by the lemma. 

The Standard Medusa. We now build a model for the mating F = f\ .u. 
/2 of the two postcritically finite quadratic maps /i, /2- We start by defining 
the standard Medusa. 

Definition 3.3. Let 0i, 62 G Z be the two rational numbers we associate to 
/1 and / 2 , as in ^2^1 Define the {61,62) standard Medusa M(6 U 6 2 ) C S 2 to 



be the union of the unit circle S 1 , the interior legs 

{pexp(2i7r2>'0i) I i<p<l,j = l,2,...} 
and the exterior legs 

{pexp(-2iir2?0 2 ) \ 1 < p < 2,j = 1, 2, . . .}. 

Defined in this way we have that z \-> 1/z maps M(#2,#i) bijectively 
to M(6 ) i,02)- The endpoints of the interior legs we denote by Xj, and the 
endpoints of the exterior legs we denote by y^, hence 

xj = 2exp(2z7r2 J ^i), j = 1, 2, . . . , and 

% = l/2exp(-2i7r2^ 2 ),j = l,2,.... 

We can think of the standard Medusa as a coupling of two standard spiders 
S>q 1 , §<9 2 , where the bodies have been cut away, then the two are glued along 
the cut. See Figure [4] for a schematic diagram of this process. 
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|z| = R |z| = R 




Figure 4. Above is a schematic of the process of mating 
1/6 with 1/7. The upper figures are the truncated spiders, 
the lower left is the Medusa on the sphere, and the lower 
right is the projection of the Medusa to the plane. 

Thurston Matings. Recall that two postcritically finite branched cover- 
ings F : § 2 — >► § 2 and g : S 2 S 2 with postcritical sets Pp and P g are called 
Thurston equivalent if there exists orientation preserving homeomorphisms 
<j) and ij) such that <j) restricted to Pp maps bijectively onto P g and o <j) 
is isotopic to the identity on S 2 rel. Pp. 

We proceed to define a branched covering p of § 2 by itself that m non- 
degenerate cases is Thurston equivalent to the mating F = f\ JL / 2 . Let 
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sIm(0i,0 2 ) be the angle doubling map rexp(z^) \-> rexp(2z0). Extend g 
smoothly to a degree two branched covering of the sphere so that: 

(1) g : D D is a degree two branched coveing with critical value at 
xi, and 

(2) g : 8 2 \B — >► 8 2 \B is a degree two branched covering with the critical 
value at y\. 

Denote by uj\ the critical point of g in B and by uj 2 the critical point of g in 
8 2 \ B. Notice that uj{ coincides with an endpoint of a leg if and only if 9{ is 
periodic under angle doubling, 9 ^ 29 mod 1. 

Notice that if we redefine g outside the unit circle to by setting it equal to 
z i y z 2 here, we obtain a map that is Thurston equivalent to f\. Similarly, if 
we instead redefine g inside the unit circle so it restricts to z \-¥ z 2 here, we 
obtain a mapping that is Thurston equivalent to f 2 . Hence it is reasonable to 
view g as our branched covering model of the mating F. Shishikura [ShiOOj 
guarantees convergence in the nondegenerate case: 

Definition 3.4. Let fi,f 2 be PCF quadratic polynomials not in complex 
conjugate limbs of M. If the two critical orbits of F = fi _u_ f 2 are disjoint, 
then fi and f 2 are called strongly mateable. 

Theorem 3.5 ([ShiOO]). i//i,/2 are strongly mateable, then g is Thurston 
equivalent to the mating F = fi JL f 2 . 

Thurston's algorithm is an iterative process that will give us a sequence 
of rational maps converging to F when F and g are Thurston equivalent. 
Using g as our model map, it works as follows. Let ao : S 2 — > C be an 
orientation preserving homeomorphism mapping ui to 0, u 2 to 00 and fixing 
1. Recursively define a n and F n as follows for n = 1,2,... Interpret o n -\ as a 
global chart defining a complex structure on S 2 . This complex structure can 
be pulled back by g. Indeed, since g is a local homeomorphism everywhere 
except at Ui, i = 1,2 we can just compose restrictions of g with o n -\. 
The complex structure defined in this way can be uniquely extended to the 
missing points uj\,U2. By the uniformization theorem S 2 equipped with the 
pullback complex structure is conformally equivalent to C. So let a n : S 2 — >► C 
be the conformal isomorphisms and normalize it so uj\ is mapped to 0, uj 2 to 
00 and 1 is fixed. By construction F n defined by the composition a n ogoa~ x 
is holomorphic. The sequence of maps constructed can be illustrated by the 
commutative diagram shown in Figure [5| 

In principle Thurston's algorithm solves our problem, the sequence of 
generated maps rational maps should converge to our mating. However the 
set of possible complex structures on S 2 is beyond actual computations, so 
we need to adapt the algorithm to allow for this. This is exactly what 
Hubbard's Medusa Algorithm does for us. 

The Medusa Algorithm. Notice that each map F n in Thurston's algo- 
rithm (in the strongly mateable described) is a degree two rational map 
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F n -{ 



s 2 — 



s 2 — 



^0 



F 2 



c 



Fi 



C 



Figure 5. A commutative diagram representing the maps 
involved in Thurston's algorithm. 



fixing 1 and having the origin and infinity as critical points. In other words, 



F n E T . By Lemma [372] we just need to know to where and oc are mapped 
to identify F n . Hence we don't need all the information contained in the 
sequence of complex structures to find F n , it is enough knowing a n -\ re- 
stricted to the standard Medusa M(#i,#2). Motivated by this we make the 
following definition. 



Definition 3.6. Set 

Mq(0 U 02) = {(T 



,02) 



a E homeo + (S 2 



and normalized}, 



where normalized here means that cr(cji) = 0,a(cj2) = oo,a(l) = 1, and 
define the Medusa space A4 as the quotient of A4o with the equivalence 
relation that identifies g\ and &2 if and only if the two maps are isotopic rel 
{xi, x 2l . . .} U {yi, 2/2, • • • , } through mappings in .Mo- 
Notice there is a natural projection tv from the complex structures on 
S 2 onto A4. Given a complex structure S we know by the uniformization 
theorem that there exists a conformal isomorphism a : (S 2 ,S) — >> C which 
we can normalize so that uo\ maps to 0, 0J2 to infinity and 1 is fixed. We let 
7r(S) equal the equivalence class of 0"|m(0i,0 2 ) m ^(^15^2)- 

One can show that there is a natural bijection between .M(0i,02) an d 
the Teichmiiller space of S 2 \ {xi, X2, . . . , yi, 2/2 5 • • • 5 1} 3 so Medusa space is 
a finite dimensional complex manifold in a natural way. 
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Mappings in Medusa space can be lifted. More precisely we have the 
following lemma. 

Lemma 3.7. Let s n -i E .Mo(0i,02) be given. Set U n — Sji — \ 

Sn-i (2/1) an( l l e ~k Fu n ,vn £ J 7 be the unique mapping as in Lemma \3J^ Then 

there is a unique mapping s n C -Mo(#i 5 #2) such that the following diagram 

commutes. 



1,02 



M(0i,0 2 r 



F Un ,v 



If s'n-1 and s n _i represent the same element in M{9i 1 02) then the lifts s' n 
and s n also represent the same element in .M(#i,#2)- 

Proof. Since the simple closed curve 7' = seperates one critical 

point and its image u n — F Un ^ Vn (0) from the other critical point 00 and 
its image v, the preimage 7 of 7' by F Un ^ Vn is a simple closed curve and 
F Un ,v n 7 —> y is a two to one covering map. Identify the fundamental group 
on S 1 with Z so that a curve having index 1 with respect to correspond 
to +1 C Z. Do similarly for 7 and 7'. Then the induced map g* : Z — » Z is 
multiplication by two. Since s n -i extends to a homeomorphism that maps 
cji to (s n -i)* : Z — Z is the identity. Finally, F Un ^ Vn maps the bounded 
component of C \ 7 onto the bounded component of C \ 7' which implies 
that (i^z n ,<y n )* : Z — >► Z is multiplication by +2. Hence (s n _i o g)* : tti(S 1 ) — >► 
7Ti(7 / ) has the same image as (i 7 ^^))* : ^1(7) — ^ follows by a 
fundamental theorem of algebraic topology that there exists a covering map 
s n : S 1 — >> 7 so that gos n -\ — s n -ioF Un ^ Vn on S 1 , and this lift is unique when 
we require that s n (l) = 1. We can extend s n to M(#i, 62) by lifting each leg 
seperately, in the way that agree with how s n is defined on the circle. In 
this way we have obtained a homeomorphism s n mapping M(#i,02) to its 
image, and we must show that s n E Mo(@i, #2)- However, since F Un ^ Vn maps 
the bounded (unbounded) part of C \ 7 to the bounded (unbounded) C \ 7' 
the image of an interior (exterior) leg is interior (exterior) , so we can extend 
s n to a orientation preserving homeomorphism of the sphere as required. 

We still need to show uniqueness of s n . For s n to be an element of 
A4o(0i,02) we must have s n (l) = 1 and that uniquely determines s n on 
S 1 . Knowing s n on the unit circle means we know to where the base point of 
the legs must lift and therefore there is only on extension to M(#i, 62) such 
that g o s n -i = s n -i o F Un ^ Vn . 

Finally suppose that and s n -i represent the same element in .M (#1, 62) 
and let s n ,s f n E .Mo(0i,02) be the two unique lifts. By assumption there 
exists an istopy connecting to s n _i, through maps in A^o(^i 5 #2)- This 
isotopy can be lifted to an isotopy connecting s n and s' n . Each map in the 
isotopy maps 1 to 1 so as before we can prove that it is an element of 

M o (e 1 ,0 2 ). □ 
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Let a starting point So E jVt(^i, ^2)5 be given. The Medusa algorithm 
consists of repeatedly applying Lemma [377] to get a sequence S n E M (61,62) 
and rational maps F Un:Vn E T for n — 1, 2, . . . . The beauty of the algorithm 
is that we produce the same sequence of rational maps that Thurston's 
algorithm produces. 

Theorem 3.8. Suppose 7r(ao) = So. Then the rational maps produced by 
Thurston's algorithm equal those produced by the Medusa algorithm, F n — 
F u n ,v n - Futhermore 7r(a n ) = S n , n = 1, 2, . . . 

Proof. Assume 7v(a n -i) — S n -i, this is case when n — 1 by assumption. 
Then s n -\ — 0"n-i|M(0i,0 2 ) * s a representative of S n -\. Now F n E T maps 
to o n -\[x\) and 00 to (T n -\(yi). So too does F Un ^ Vn . Hence, by Lemma 3.2 
F n = F UniVn . We have that <7 n | M(g g a) E Mo(6 Xl 2 ), is a lift of s n _i. So by 



the uniqueness part of Lemma 3.7 7r(cr n ) = 5 n . The theorem now follows by 



induction. □ 



Now we can justify the Medusa algorithm, by combining Theorems 2.1 



3.5, and 3.8 



Theorem 3.9. If f\ and f 2 are strongly mateable, then the Medusa algo- 
rithm converges to the mating F = f\ _u_ f 2 . 

In practice, the algorithm seems to converge without assuming the maps 
are strongly mateable. Thus we expect that a stronger theorem holds; 
namely, it should be the case that anytime f\ and f 2 are PCF quadratic 
polynomials in complex conjugate limbs of M, the Medusa algorithm should 
converge to the mating. The case not covered by Thurston's theorem is when 
two polynomials that are not in complex conjugate limbs have a mating with 
only one critical orbit. In this case naively running the Medusa algorithm 
produces a sequence of Medusas which does not converge (rather tends to 
the boundary of the Teichmuller space), but the obstruction points (the 
critical orbits becoming identified) are all pushed together upon iteration of 
the algorithm, hence the sequence of rational maps seems to converge to the 
mating. To prove this stronger result one could investigate how the maps 
in the Medusa algorithm are converging as the boundary of the Medusa 
space is approached. We expect the techniques of Nikita Selinger's PhD 
thesis [SellOj on convergence at the boundary of Teichmuller space could be 
adapted to solve this question, and leave this future result to the interested 
reader. 

3.2. The Implementation. The point of the Medusa algorithm is that it 
lends itself to implementation as a computer program. The implementation 
is an adoption of the implementation of the spider algorithm to the more 
general setting of quadratic rational maps. 

To initiate the program, the user inputs two rational angles 6\ 1 6 2 . The 
implementation defines an initial Medusa so'- M.(6i,6 2 ) C, say close to 
the identity. 
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To describe our matings, we define a chart on T by letting R a ^ : z \-> 
^z 2 X\i-b) 1 ( a '^) ^ \ {( z i z ) I z ^ Q- I n this wa y we parametrize all the 
maps in J 7 . Supposing that F E T maps to ^ and oc to v, we let a = V ^Z^ 
and 6 = g=i. Then i? a , 6 = F = F u , v . 

We represent a mapping 5 : M(#i, #2) — > C by several lists of points in C. 
One list represent the image of the unit circle, and the other lists represent 
the images of the legs. Also we always let the list of points representing the 
image of the unit circle start with the point 1. 

We adopt the convention that two consecutive points in the image of 
the unit circle or in a leg is connected by an arc of circle. For the points 
on the image of the circle or on the interior legs the circle chosen is that 
through and the arc of circle chosen is the one connecting the two 

points and omitting s(yi). For consecutive points on the exterior legs adopt 
the convention that they are connected by the arc of the circle through the 
points and s(xi). The arc is the one that connects the two points and omits 
s(xi). 

Clearly, with the information contained in the lists of points and the 
convention just mentioned we can reconstruct, not s, but the isotopy class 
of s. 

An iteration consists of finding the class of the pullback of s n -i (as in 



Lemma 3.7). As in the implementation of the spider algorithm we break 
the process down into three steps: a pullback step, a rectifying step and a 
pruning step. 

Pullback. Given s n -i as lists of points as described we first find F Un ^ Vn = 
R an ,b n - This corresponds to solving 

(2) ^— ^ = u n = s n -i(xi) and ^ = v n = s n _i(yi). 



In other words 



(u n -l)v u n -l 
and b n 



u—v u—v 
Notice that R an ,b n is the composition of a Mobius transformation with z 1— >> 
z 2 . Hence, pulling back a point consists of first pulling it back by a Mobius 
transformation M n and then by the square. The question that needs to 
resolved is, what branch of the squareroot do we need to choose. 

First we pullback the points corresponding to the image of the unit circle. 
Suppose that we have pulled back a point z^ and obtained the point Wk 
and want to pullback the next point in the list Zk+i- Pulling back first by 
the Mobius transformation we get that the circle through and v n 

becomes a circle through M~ 1 (z/ C ), M~ 1 (z/ c+ i) and oc i.e. a line. Since the 
arc of circle connecting the two points was chosen to be the one that did not 
contain v n the pullback of the arc of circle by the Mobius transformation 
becomes simply a line segment between M~ 1 (z^), M~ 1 (z/ c+ i). The preimage 
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of a line by the square is a hyperbola, the two branches of which are con- 
tained in opposite quarter planes. Hence knowing one preimage Wk, we need 
to choose the square root so that and + 1 lies in the same halfplane. 

So the pullback the points corresponding to the circle we construct to 
lists, A, B. The first element of A is 1 and the first element of B is — 1, i.e. 
the two preimages of 1 by R an ^ n . This was the first step. Next we iterate 
through the remaing points in the list. The fc'th step consists in finding the 
two preimages of call them and w' k . If the last inserted point in the 
list A lies in the same quarter plane as Wk then we insert Wk in A, and w k 
in B. Otherwise we insert Wk in B and w k in the A. It is easy to verify that 
the points in the list A are images of the points on the unit circle of angles 
in the interval < 9 < 7r, whereas the points in B correspond to angles 9 
with TV < 9 < 2tv. Having pulled back all the points we can concanate the 
two list so we get one list (starting with the point 1) representing the image 
of the circle by s n . Notice that this list contains twice the points of the one 
we have just pulled back. 

Next we pull back the interior legs. The leg corresponding to angle 9 is the 
preimage of the leg corresponding to angle 29. If < 9 < tt the point in the 
list A list that is the preimage of the anchor point of the leg of angle 29 will 
be the anchor point of the new 9 leg, otherwise it will be the corresponding 
point in the list B. Hence we have already computed (and can localize) the 
pull-back of the first point in the leg. Hence, as before we can pull back the 
rest of the leg, we need to chose the square root so the consecutive points 
lies in the same halfplane. 

Pulling back the outer legs is essentially the same, except that now pulling 
back by M n two consecutive defines a arc of circle, where the circle goes 
through 0. However since z \-> z 2 commutes with z \-> l/z we can write the 
squre as the composition of l/z, z 2 and then l/z again. Hence pulling back 
by M n and then making the change of coordinates w = l/z we are back in 
the same situation as the one we were facing when pulling back the interior 
legs. 

In this way we obtain a list of point representing the map s n . However, 
the points are now connected by arcs of hyperbolas and not arcs of circles. 
The next step, rectifying, remedy this situation. 

Rectifying. Perhaps a better word for the second part of an iteration would 
be circlyfying. We want to bring us back to the starting position where 
consecutive points in the lists are connected by arcs of circles. This is the 
most delicate part of the implementation. What we want to do is replace the 
arcs of hyperbolas with arcs of appropiate circles without changes the isotopy 
class of the corresponding element in M.{9\ 1 92)- So given two consecutive 
points zi, Z2 we want to see if there is a homotopy from an arc of hyperbola 
to an arc of circle so that the intermediate curves does not cross any of 
the distinguished points s n (xi), s n (x2), s n (2/2)5 •• •}• It is rather 

tedious so we will only outline how it is done. The circle and the hyperbola 
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are two (real) quadratic curves and we first find their intersection. This 
can be boiled down to finding the roots of a degree 4 equation in one real 
variable. However, since we know that z\ and Z2 lies on both curves, we can 
do a division of polynomial and the remaing points (if any) can be found 
by solving a quadratic equation. The most difficult case when the branch of 
hyperbola containg z\ and Z2 intersect the circle in four points. Then the 
union of the circle and the branch of hyperbola cuts the plane into six parts. 
By elementary geometric reasoning, one can find exactly to which of the six 
parts a given point belongs, and this knowledge is enough to decide if the 
homotopy exists. 

If the homotopy exists then we can move on, but if it doesn't we need to do 
something. What we do is to subdivide the arc of hyperbola in two halves, 
ziX and £, Z2 and recursively rectify each half. In case we are not dealing 
with a leg terminating at a distinguished point, then by compactness the 
distinguished points are a definite distance away from the arc of hyperbola 
between zi,^2- Given any e > 0, any fine enough subdivision of the arc 
of hyperbola, zi : (j, (2, • • • , (k, z 2 will satisfy, that if we replace the parts of 
hyperbolas with arcs of circles we will stay with a spherical e neighborhood of 
the original arc of hyperbola. Hence, we are able to rectify after adding only 
a finite number of points. In the case that the arc of hyperbola terminates 
in a distinguished point Z2 then we are dealing with the image of a leg. It is 
not difficult to see that we do not change the isotopy class of s n by allowing 
the homotopy to cross Z2. In practice, this means that when rectifying a leg, 
we do not consider the endpoint of the leg a distinguished point, and we are 
sure that we can rectify adding only a finite number of points. 

Pruning. After pulling back and rectifying, we have new lists of point rep- 
resenting 5 n , but the number of points representing the image of the unit 
circle has at least doubled. This means that unless we do something we will 
run out of memory in a finite number of iterations. 

What we do is pruning which amounts to checking each point Z2 that is not 
the attachment point or terminal point of the leg whether it can be removed 
without changing the isotopy class of the represented map. In practice this 
means checking whether two arcs of circles, one through z\ and Z2 the other 
through Z2 and £3, can be replaced by an arc of circle going from z\ to z% 
without changing isotopy class. Using a Mobius transformation to change 
coordinates the question becomes whether a line segment (^1,^2) and a 
line segment (u>2, ^3) can be homotopied to a line segment (u>i, W3) without 
crossing distinguished points, a question that can be easily answered. 

Drawing the Julia set. In addition to producing a sequence of maps 
Ra n ,b n converging to the mating, the Medusa algorithm can be used to draw 
successive approximations to the Julia set of the mating. At the beginning of 
the program, a "painted" sphere /Co is created, with each point in the upper 
hemisphere painted black, and each point in the lower hemisphere painted 
white (or clear). At each iteration of the algorithm, given parameters a m , b m 
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Figure 6. Each of the three columns above shows Maple 
output of the actual Medusas used in the iteration of the 
Medusa algorithm for the mating of 1/7 with 1/3 (rabbit 
mate basilica). In each column, the top figure is the Medusa 
on the sphere, the lower figure is the Medusa projected onto 
the plane. Leftmost is the initial Medusa, central is after 2 
steps, rightmost is after 20 steps. 



and a painted sphere K m -\ (i.e., a sphere with each point marked one of 
black or white), the program computes the pull back of K m -\ by it!" 1 & , to 
create JC m . 

When the sequences (a m ,6 m ) converge, then R am ,b m converges to R a ^ = 
fi _u_ / 2 , and K m converges to /C, with white or clear marking the Julia set 
of /i, and black the Julia set of /2. 

For example, let C1/4 be the parameter which is the landing point in the 
Mandelbrot set of the external ray of angle 1/4 (ci/4 « —0.228 + 1.115i). 
This is a tip point on the rabbit bulb. The mating of z 2 + c^u with itself 
exists, and is studied in detail in [Mil04j. In this case the Julia set of the 
mating is the entire sphere, so the approximations K n drawn by Medusa are 
particularly interesting. Figure [7] shows approximations Kq, /C10, and /C14 
for this mating. Also see §4.3| for other similar examples. 

The full source code for Medusa is available for download at |Dyn| . There 
are still a few bugs, most notably: when mating with a p/q where q is even, 
the algorithm will converge properly for a few steps, then start diverging. 
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Figure 7. Upper Left: the Julia set of f(z) = z 2 + cy^ 
where C1/4 is the landing point in the Mandelbrot set of the 
external ray of angle 1/4 (so c^u ~ —0.228 + 1.115i), shown 
with critical orbit rays 1/4, 1/2,0. Clockwise around (upper 
right, lower right, lower left): approximations on the sphere 
Kq, /C10, /C14, respectively, to the Julia set of f(z) = z 2 + C1/4 
mated with itself. 

4. Examples 

In this section we discuss several types of matings with different proper- 
ties. For simplicity, we will refer to a PCF quadratic polynomial simply by 
its rational angle 9 C — p/q, or sometimes f p / q . 

4.1. Simple examples. We explain our first example of an image of a 
mating produced by the Medusa algorithm in detail. We will mate the two 
quadratic polynomials shown in Figure [8} fa will be the rabbit, 1/7, and fa 
will be the basilica, 1/3. 

Let F — fa AL f 2 = 1/7 JL 1/3. The rightmost sphere in Figure [8] illustrates 
the Julia set of the mating F. Due to our normalization (Equation ??), 
the critical point of f\ is always at z — in the sphere, shown as the 
south pole, and the critical point of fa is sent to z — 00 in the sphere, 
shown as the north pole. The portion of the filled Julia set of the mating F 
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Figure 8. From left to right: The Julia set of the rabbit, 
critical angle 1/7, then the Julia set of the basilica, critical 
angle 1/3, both shown with both sets of critical orbit rays 
(1/7,2/7,4/7,1/3,2/3) for comparison; finally, the mating 
1/7 mate 1/3 on the sphere, with 1/3 in black, and 1/7 clear. 



which corresponds to J{f\) (the rabbit) is shown in clear, and "centered" 
about the north pole. The portion corresponding to J(f2) (the basilica) is 
shown in black on the front half of the sphere, and grey on the back half 
(to indicate that to see this, you are looking through J(/i)). However, due 
to the symmetry of the Julia sets of quadratic polynomials, this image is 
invariant under 180 degree rotation about the vertical axis, hence the grey 
image in the back does not convey new information. Also, the fixed point 
z = 1 (corresponding to the /3-fixed points of /i,/2), is in the dead center 
of the image, in the front. Note reversing the order of mating, drawing the 
image of 1/3 JL 1/7, would have the effect of a 180 degree rotation about the 
central horizontal axis (from z = 1 to z = — 1), and flipping the colors. 

Self-mating. The limb of the mandelbrot set enclosed by rays of angle 
1/3, 2/3 (see Figure [3]) is the only limb which is its own complex conjugate. 
As such, any PCF quadratic polynomial which is not in that limb can be 
mated with itself. Such a mating clearly has extra symmetries. The leftmost 
image in Figure [9] is the rabbit 1/7 mated with itself. We discuss self matings 



more in £4.4 



Tuning. One simple way to make a mating more complicated is by tuning 
one of the quadratic polynomials. The result shows up as you would expect. 
In figure [9j compare the rabbit mate rabbit on the left with the right figure, 
in which the clear rabbit has been tuned with a basilica. We explore further 



expectations (and surprises) concerning tunings in £4.5 



4.2. Shared Matings. One of the intriguing observations in the study of 
matings is that it can happen that two distinct pairs of PCF quadratic 
polynomials give rise to the same mating F. If f\ JL / 2 = F = fs JL /4, and 
fi 7^ H or /2 / /4, then we call F a shared mating. 
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Figure 9. Left: 1/7 mate 1/7, rabbit mate rabbit. Right: 
1/7 mate 10/63, i.e., replace each disk in the leftmost clear 
rabbit with a basillica. 



The simplest kind of shared mating is when f\ _ll / 2 = 11 fi- For 
example, the left side of Figure [10] illustrates such a shared mating of the 
rabbit (1/7) and aeroplane (3/7). Of course, taking a shared mating and 
performing the same tuning on each quadratic polynomial will produce an- 
other shared mating, for example as on the right side of Figure [10} 

Wittner ( | Wit 88] ) studied this, and related, shared matings. 

4.3. Space-filling curves and Lattes mappings. A very different exam- 
ple of a shared mating, discussed in detail in |Mil04] . is a Lattes map which 
can be realized as a mating in four distinct ways: 

1/6 JL 5/14 ^ 3/14 jl 3/14 ^ 3/14 jl 1/2 ^ 5/6 JL 1/2. 

It is not known whether there is a bound on the number of ways in 
which a quadratic rational map can be realized as a mating. The quadratic 
polynomials involved above are: /i/q(z) = z 2 + z, a tip point on the rabbit 
limb; f^/e(z) = z 2 — i, the complex conjugate of fi/Q] /5/14, a tip point of the 
bulb on the basilica bulb corresponding to the rabbit; and f 1/2(2) = z 2 — 2, 
the real tip point of the basilica limb (the leftmost point in the mandelbrot 
set). The Julia set for each of 1/6, 5/14, 3/14, 1/2 is a dendrite, hence has 
empty interior. For example, the Julia set of fi/± is a dendrite, shown in 
Figure [7| Below is a characterization of when this occurs. 

Fact 4.1. Suppose P c is a PCF quadratic polynomial. Let 9 C = p/q be a 
reduced fraction. TFAE: 

(1) K c has empty interior; 

(2) q is even; 

(3) 9 C is strictly pre-periodic under angle doubling. 



Thus the mating of any two quadratic polynomials satisfying Fact 4T 
(including the shared mating above) has Julia set the entire Rieman sphere. 
You can visualize such a mating as a space-filling curve on the sphere (each 
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Figure 10. Upper left: the rabbit, 1/7; Upper right: the 
aeroplane, 3/7, both shown with both sets of critical or- 
bit rays (1/7,2/7,4/7,3/7,6/7,5/7). Lower left: the shared 
mating the rabbit mate the aeroplane, 1/7 mate 3/7, equiva- 
lently, the aeroplane mate the rabbit. Lower right: basilicas 
in the rabbit mate basilicas in the aeroplane, 10/63 mate 
28/63. 

of the empty interior Julia sets is a curve which is pulled into becoming a 
space-filling curve). Further, since the Julia set of fi/ 2 is a line segment, 
any mating of the form p/q _u_ 1/2 where q is even will create a space- filling 
Peano curve. 

Since the Julia set is the entire Riemann sphere, we cannot very well study 
such matings by drawing their Julia sets. The harmonic measure supported 
on the Julia set is an object which deserves further study. One could hope to 
learn something by examining the approximations to the Julia set drawn by 
the program Medusa in the steps of the algorithm converging to the mating. 
See Figure [TTJ 

4.4. Self Matings. Carston Peterson has observed that if / is any PCF 
quadratic polynomial which is not in the 1/2-limb of the Mandelbrot set 
(i.e., not in the unique limb which is its own complex conjugate), then the 
following two rational maps are topologically conjugate: 
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Figure 11. Each of the four images above illustrates a 
Medusa approximation K\2 to the same shared Lattes mat- 
ing. Upper left: 1/6 mate 5/14. Upper right: 3/14 mate 
3/14. Lower left: 1/2 mate 3/14. Lower right: 1/2 mate 
5/6. (Note the two lower figures are mated in reverse order 
from the shared mating. Just rotate the picture 180 degrees 
and exchange the colors to see the correct image). 



(1) start with / _u_ /, then mod out by the obvious symmetry, and 

(2) / ^ A/2, where f 1/2 (z) = z 2 -2. 

This is because for f 1 / 2 ^ the Julia set is a line segment, [—2, 2], and every 
external ray of angle 9 has the same landing point as the ray of angle 1 — 9 
(the ray is horizontal and lands at 2, the ray of angle 1/2 is horizontal and 
lands at —2). 

For example, shown in Figure 12 is the Julia set of A/5, together with the 
Julia sets of both the self mating of /1/5, and the mating of 1/5 with 1/2. 
Since the Julia set of 1/2 is simply a line segment, note in the figure how 
this simple segment is twisted to fill up all of the black. 

4.5. Sequences of matings, and their limits. One question about mat- 
ings which has yielded an interested study is: If f\ and f 2 are quadratic 
polynomials not in complex conjugate limbs, which are not PCF, when does 
a mating exist (assuming connected Julia sets)? If A and f 2 are hyper- 
bolic, thus stable perturbations of hyperbolic PCF polynomials 51,52? each 
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Figure 12. Upper Left: The Julia set of f ± / 5l which is the 
center of largest baby Mandelbrot set off of the rabbit bulb, 
shown with critical orbit rays 1/5,2/5,4/5. Upper Right: 
1/5 mate 1/5. Lower Right: 1/2 mate 1/5, i.e., mod out 
the upper figure by the obvious symmetry. Lower Left: an 
approximation /Ci6, to 1/2 mate 1/5. The black is 1/2, so 
shows the simple line twisting to fill up the alloted space. 



with a super attracting periodic cycle, the mating exists as a deformation 
of the mating of 51,52- Several papers have appeared constructing matings 
between particular non- hyperbolic polynomials (see Haissinksy and Tan Lei 
|HL04| . Luo [Luo95| . Yampolsky, Zakeri jYZOlj .) However, Epstein [E^s] 
has shown that mating does not extend continuously to the boundary of the 
hyperbolic component (in fact, the set of points in dM x dM where there is 
no continouous extension is dense). Epstein's theorem is that an obstruction 
to continuously extending this map to a mating between the two root points 
of the hyperbolic components occurs whenever in the mating 51 JL 52, the 
immediate basins of the superattracting cycles of 51, 52 touch along a distin- 
guished repelling cycle (excluding gi{z) — z 2 ). For example, this occurs in 
the mating of the rabbit and the aeroplane, Figure [T0| That this is a shared 
mating is an additional coincidence, not needed for Epstein's theorem. 
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Figure 13. Medusa images of the Julia sets of the follow- 
ing matings: Left: 1/511 mate 1/511. Center: 1/511 mate 
255/511. Right: 169/511 mate 169/511. I should elaborate... 



We can use Medusa to see a different type of example of why mating as a 
map from MxMto the space of quadratic rational maps is not continuous. 
We examine a few convergent sequences of quadratic polynomials, m , uj m —> 
0, cj, as m —> oo, such that the mating 9 m JL uj m exists for every m. but 9 iL 
uj either does not exist, or is not the limit of 9 m _u_ uj m . 

Below are some simple examples of sequences with no limit, or the wrong 
limits. 

(1) First consider 9 m — uo m — ? 9 — u — 0. Note corresponds 
to z \-> z 2 , so 9 ii. cj = ii. is just z z 2 , with Julia set the circle. 
However, Medusa output suggests that the Julia set of 9 m _u_ uj m is 
much more complicated than the unit disk. The leftmost image in 



Figure 13 shows the Julia set of 1/255 ii- 1/255 (recall Figure [9] shows 
the first element of the sequence, 1/7 JL 1/7). 
(2) A similar example is given by 9 m — 2m 1 _ 1 -> u rn — ^mj^ 1 ? so = 
and uj — 1/2. Note JL 1/2 is just 1/2, i.e., z \-> z 2 - 2, with Julia 
set [—2,2]. As in the previous example, Medusa output shows 9 m 



JL (jj m is quite complicated. The center of Figure [13] shows 1/511 JL 
255/511 (and Figure [l0| shows the first element of the sequence, the 
rabbit mate the aeroplane). 

(3) Finally, we examine 9 m — uj m — ^ 2 2^+3 2 /i^~ 1 (i- e -> ^ ne sequence 
9/31,41/127, 169/511, . . ., of angles of the upper ray landing at the 
root point of the bulbs proceeding from the rabbit to the basilica), 
hence 9 = uj = 1/3. Since fi/s(z) = z 2 — 1 is the basilica, it is 
its own complex conjugate, and its self mating does not exist. The 



rightmost image in Figure 13 is the Julia set of 169/511 mated with 
itself. 
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